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Abstract 

In the data of the ether-drift experiments there might be sizeable fluctuations superposed 
on the smooth sinusoidal modulations due to the Earth's rotation and orbital revolution. 
These fluctuations might reflect the stochastic nature of the underlying " quantum ether" and 
produce vanishing averages for all vectorial quantities extracted from a naive Fourier analysis 
of the data. By comparing the typical stability limits of the individual optical resonators with 
the amplitude of their relative frequency shift, the presently observed signal, rather than being 
spurious experimental noise, might also express fundamental properties of a physical vacuum 
similar to a superfluid in a turbulent state of motion. In this sense, the situation might be 
similar to the discovery of the CMBR that was first interpreted as mere instrumental noise. 



1. Introduction 



In principle, the physical vacuum might behave as a medium with a non-trivial refractive 
index. In this case, the speed of light in the vacuum c~ = -jt-^ — would not coincide ex- 

' ^Vvacuum 

actly with the parameter c entering Lorentz transformations and therefore could propagate 
isotropically in only one frame X. In a moving frame S', there would be an anisotropy of the 
two-way speed of light c 7 (#) 

c c z 

with p] 

\B\ ~ 3(AA vacuum - 1) (2) 

V being the speed of S' with respect to E. This anisotropy could then be measured through 
the beat frequency Au of two orthogonal cavity-stabilized lasers. 

In this context, the search for time modulations of the signal induced by the Earth's 
rotation (and its orbital revolution) has always represented a crucial ingredient for the analysis 
of the data. For instance, let us consider the relative frequency shift of two optical resonators 
for the experiment of Ref.|2j 

= S(t) sin 2uj Iot t + C(t) cos 2uj Iot t (3) 

where a; ro t is the rotation frequency of one resonator with respect to the other which is 
kept fixed in the laboratory and oriented north-south. If one assumes that, for short-time 
observations of 2-3 days, the time dependence of a hypothetical physical signal can only be due 
to (the variations of the projection of V in the interferometer's plane caused by) the Earth's 
rotation, S(t) and C(t) admit the simplest Fourier expansion (r = co s idt is the sidereal time 
of the observation in degrees) [2] 

S(t) = So + S s i sinr + S c i cos r + S S 2 sin(2r) + S C 2 cos(2r) (4) 

C(t) = C + C s i sin r + C cl cos r + C s2 sin(2r) + C c2 cos(2r) (5) 

with time-independent and coefficients. Therefore, by accepting this theoretical frame- 
work, it becomes natural to average the various Ck and Sk over any 2-3 day period. By 
further performing inter-session averages over many short-period experimental sessions (see 
Fig. 2 of Ref.[3]) the general conclusion [3] is that, although the typical instantaneous signal 
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is 0(1O~ 15 ) or larger, the average Ck and S& coefficients are at the level 0(1CT 17 ) and, with 
them, the derived parameters entering the SME [5] and RMS [5] models. 

However, as we shall argue in Sect. 2, in an alternative picture of the vacuum, one could 
consider different types of ether-drift. In this case, the same experimental data could admit 
a completely different interpretation. This possibility will be illustrated in Sects. 3 and 4 
where, by comparing with the typical stability limits of the individual optical resonators, we 
shall argue that the presently observed signal is unlike to represent just spurious instrumental 
noise. For this reason, in Sect. 5, we shall discuss a possible theoretical framework that can 
explain the present data. In this sense, as discussed in the final Sect. 6, the situation with the 
ether-drift experiments might resemble the discovery of the CMBR that, at the beginning, 
was also interpreted as mere instrumental noise. 

2. The vacuum as a turbulent superfluid 

By accepting the idea that there might be a preferred reference frame, which is the modern 
denomination of the old ether, one could also change perspective and, before assuming any 
definite theoretical model, first ask: if light were really propagating in a physical medium, 
an ether, and not in a trivial empty vacuum, how should the motion of (or in) this medium 
be described ? Namely, could this relative motion exhibit variations that are not only due to 
known effects as the Earth's rotation and orbital revolution ? 

Without fully understanding the nature of that substratum that we call physical vacuum, 
it is not possible to make definite predictions. Still, according to the standard model of 
electroweak and strong interactions, this physical vacuum could be considered as a superfluid 
medium [6] filled by those particle condensates that play a crucial role for many fundamental 
phenomena such as mass generation and quark confinement. By further considering the 
idea of a non-zero vacuum energy, this physical substratum could also represent a preferred 
reference frame [7J. In this picture, the standard assumption of smooth sinusoidal variations 
of the signal, associated with the Earth's rotation and its orbital revolution, corresponds to 
describe the superfluid motion in terms of simple regular flows. 

However, visualization techniques that record the flow of superfluid helium show [8] the 
formation of macroscopic turbulent structures with a velocity field that fluctuates randomly 
around some average value. The same idea of turbulence, that dates back to the original Lord 
Kelvin's attempts to generate transverse waves in the ether [9], is also suggested by the formal 
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equivalence that one can establish [lOl E] between the propagation of small disturbances in 
a turbulent fluid and the propagation of electromagnetic waves as described by Maxwell 
equations. 

To exploit the possible implications for the ether-drift experiments, let us first recall the 
general aspects of a turbulent flow. This is characterized by extremely irregular variations of 
the velocity, with time at each point and between different points at the same instant, due to 
the formation of eddies [12J. For this reason, the velocity continually fluctuates about some 
mean value and the amplitude of these variations is not small in comparison with the mean 
velocity itself. The time dependence of a typical turbulent velocity field can be expressed as 

n 

v(x,y,z,t)= %iP2.4>™(^y> 2 ) ex P( _i X^'^') ( 6 ) 

PiP2--Pn j = l 

where the quantities = ujt + /3j vary with time according to fundamental frequencies ujj 
and depend on some initial phases (5j. As the Reynolds number 7Z increases, the total number 
n of ujj and (5j increases. In the 1Z —* oo limit, their number diverges so that the theory of 
such a turbulent flow must be a statistical theory. 

As discussed by Feynman [13], the idea that there might be a rich structure of turbulent 
eddies for very large Reynolds numbers is, at first sight, quite surprising. After all, starting 
from the Navier-Stokes equation, the eddies are generated, in the equation for the vorticity 
ft = V x v, by the viscous term which formally vanishes in the 1Z — > oo limit. However, there 
is a subtlety. The right-hand term of the equation 

^- + Vx(f!xv) = -An (7) 

has ^- times a second derivative. This is the higher derivative term in the equation. Thus, 
although ^ becomes smaller and smaller, there can be rapid variations of ft, producing larger 
and larger Af2, that compensate for the small coefficient. For this reason, the solutions of 
the viscous equation do not approach the solutions of the equation 

— + V x (n x v) = (8) 

when ^ — > 0. 

Now, due to the presumably vanishingly small viscosity of a superfluid ether, the relevant 
Reynolds numbers are likely infinitely large in most regimes and we might be faced precisely 
with such limit of the theory where the physical vacuum behaves as a stochastic medium. 
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In this case random fluctuations of the physical signal, superposed on the smooth sinusoidal 
behaviour associated with the Earth's rotation (and orbital revolution), would produce devi- 
ations of S(t) and C(t) from the simple structure in Eqs.flU) and © and an effective temporal 
dependence of the fitted = Ck(t) and Sk = Sk(t) vectorial coefficients. In this situation 
one could easily get, due to phase interference, vanishing average values (Ck) = (Sk) = 0. 

Nevertheless, as it happens with random fluctuations, the average amplitude of the signal 
could still be preserved. Namely, by extracting from the data the positive-definite combination 

A(t) = VS 2 (*) + C 2 (t) (9) 

a definite non-zero (A) might well coexist with (Ck) = (Sk) = 0. 

3. Noise or stochastic turbulence ? 

To provide some evidence that indeed we might be faced with this type of situation, let us 
consider the experimental apparatus of Ref.|14j where the two optical cavities were obtained 
from the same monolithic block of ULE. Due to sophisticated electronics and temperature 
controls, the stability limits for the individual optical cavities are extremely high. More 
precisely, the effect of residual amplitude modulation is below 0.02 Hz for both cavities and 
thus about 7 • 10 -17 for a laser frequency 2.82 • 10 14 Hz. For the non-rotating set up, and with 
power stabilization, the laser power fluctuations are below 3 • 10~ 17 for cavity 1 and below 
1 • 10~ 17 for cavity 2. Finally, for the non-rotating set up, also the tilt instabilities are below 
1 • 10~ 16 , or 0.03 Hz, for both cavities. Thus, by adding all effects, one deduces a stability of 
about ±0.05 Hz for the individual cavities. This is of the same order of the mean frequency 
shift between the two resonators, say (Au) ~ ±0.06 Hz, when averaging the signal over long 
temporal sequences. 

However, the instantaneous Au is much larger, say ±1 Hz, and so far has been interpreted 
as spurious instrumental noise. To check this interpretation, we observe that, in the absence 
of any light anisotropy, the noise in the beat frequency should be comparable to the noise of 
the individual resonators. Instead, for the same non-rotating set up, the plateau of the Allan 
variance for the beat signal was found 10 times bigger, namely 1.9 -10~ 15 with a corresponding 
typical shift of about 0.56 Hz (see Fig. 8 of Ref. |14j ) . Also the slopes are different in the two 
cases, indicating noises of different nature. The authors tend to interpret this as cavity 
thermal noise and refer to [15] . However, if the theoretical estimate of Ref.[15j applies, the 
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relevant effect should be considerably smaller, about 4 • 10~ 16 or 0.13 Hz (see the first entry 
in Table I of [IS])- I n anv case, one can also compare with experiments performed in the 
cryogenic regime. If this typical 0(1O~ 15 ) beat signal reflects the stochastic nature of an 
underlying quantum ether, it should be found in these different experiments as well. 



4. An alternative analysis of the data 

Motivated by the previous arguments, we have explored the idea that the observed beat signal 
could be due to some form of turbulent ether flow. In principle, in this perspective, one should 
abandon the previous type of analysis based on a fixed preferred reference frame and extract 
the average amplitude from the instantaneous data before any averaging procedure. However, 
these instantaneous data are not available and thus, to extract (A), that was never reported 
by the experimental groups, one can only use the C\. and Sk coefficients averaged within each 
2-3 day session. Due to the expected negative phase interference, the resulting (A) should 
represent a lower limit for its true experimental value. 
To extract (A), one should first re- write Eq.([3]) as 

= A(t) cos{2oj TOt t - 20 o (*)) ( 10 ) 

with 

C{t) = A(t) cos 20 (t) S(t) = A(t) sin 29 (t) (11) 

9o(t) representing the instantaneous direction of the ether-drift effect in the plane of the 
interferometer. In this plane, the projection of the full V is specified by its magnitude 
v = v(t) and by its direction 9q = &o(t) (counted by convention from North through East so 
that North is 9q = and East is 9q = 7r /2). If one assumes Eqs.(jl]) and ([5]), then v(t) and 
6o(t) can be obtained from the relations [16\ [T7] 

cos z{t) = sin 7 sin + cos 7 cos <f> cos(r — a) (12) 

sin z(t) cos 0o(t) = sin 7 cos (ft — cos 7 sin (f> cos(r — a) (13) 

sin z(t) sin 9o(t) = cos 7 sin(r — a) (14) 

v(t) = Vsinz(t), (15) 

where a and 7 are respectively the right ascension and angular declination of V. Further, <p 
is the latitude of the laboratory and z = z(t) is the zenithal distance of V. Namely, z = 
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corresponds to a V which is perpendicular to the plane of the interferometer and z = vr/2 to 
a V that lies entirely in that plane. From the above relations, by using the 0(v 2 /c 2 ) relation 
A(t) ~ 4^' the other two amplitudes S(t) = A(t)sm29 (t) and C(t) = A(t) cos 26 (t) can 
be obtained up to an overall proportionality constant. By using the expressions for S(t) and 
C(t) reported in Table I of Ref. [2 J (in the RMS formalism [5]), this proportionality constant 
turns out to be \\B\ so that we finally find 

Att) = \\B\^ (16) 
2 c z 

where B is the anisotropy parameter entering the two-way speed of light Eq. (fTJ) . It is a simple 
exercise to check that, by using Eqs. (fTTj) . Eqs. (fT2"j) - (fTl)j) and finally replacing % = 90° — (f>, one 
re-obtains the expansions for C(t) and S(t) reported in Table I of Ref. [2]. 

We can then replace Eq. (fT5|) into Eq. (fl6|) and, by adopting a notation of the type in 
Eqs.©-©, express the Fourier expansion of A(t) as 

A(t) = A + Ai sin r + A 2 cos r + A 3 sin(2r) + A A cos(2r) (17) 

where (the daily averaging of any quantity is here denoted by (..)), 

(A) = Aq = ^\B\^p- = \\B\^ (l - sin 2 7 cos 2 x- ^ cos 2 7 sin 2 X ) (18) 



IV 2 IV 2 

A\ = — — ^- sin 27 sin a sin 2% A2 = — — \B\—=- sin 2 7 cos a sin 2\ (19) 

4 c z * 4 cr 

l^,^ 2 2 2 - l.~.^ 2 



^3 = — -rl-Bl— 7T cos 7sin2asin x A4 = —-\B\—^- cos 7cos2asin x (20) 

A c z 4 cr 

To obtain from the and S^, we observe that by using Eq, (|17|) one obtains 

(A 2 (t))=A 2 + ±(A 2 +A 2 + A 2 + A 2 ) (21) 
On the other hand, by using Eqs.flU), © and (fTTI) . one also obtains 

( A 2 {t) ) = { C 2 (t) + S 2 (t) ) = Cl + S 2 Q + Q 2 (22) 

where 

Q = \H(C 2 U + S 2 U + C 2 2 + S 2 2 ) (23) 

and 
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S u = yj + S 2 2 = \J ^ + S c2 

Therefore, one can combine the two relations and get 



(25) 



Al{l + r) = Cl + Sl + Q 2 



(26) 



where 




(27) 



< r < 0.4 



(28) 



for the latitude of the laboratories in Berlin [2] and Diisseldorf [IB] in the full range < | —/ 1 < 
7r/2. We can thus define an average amplitude, say Aq, which is determined in terms of Q 
alone as 

A = - = ~ (0.92 ± 0.08)Q (29) 
yl + r 

where the uncertainty takes into account the numerical range of r in Eq. (|28|) . This quantity 
provides, in any case, a lower bound for the true experimental Aq since 



At the same time Q is determined only by the C s \, C c \ r .. and their S-counterparts. According 
to the authors of Refs.[2J[3], these coefficients are much less affected by spurious effects, as 
compared to Cq and Sq, and so will be our amplitude Aq. 

We have computed the Q values for the 27 short-period experimental sessions of Ref.[3]. 
Their values are reported in Table I. Although more recent data were announced in Ref.jlj, 
only the results of fits within the SME model, and not the data themselves, are explicitly 
reported. In any case, the data of Ref.[3] represent, within their statistics, a sufficient basis 
to deduce that a rather stable pattern is obtained. This is due to the rotational invariant 
character of Q in the 8-th dimensional space of the C s i, C c i,...S S 2, S C 2 so that variations of the 
individual coefficients tend to compensate. By taking an average of these 27 determinations 
one finds a mean value 




(30) 



Q = (13.0 ±0.7 ±3.8) • 10 



16 



Ref.[3] 



(31) 
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where the former error is purely statistical and the latter represents an estimate of the sys- 
tematical effects. Such systematic uncertainty, deduced from the data of Ref.[3], is also 
completely consistent with the comparable stability limit 4- 1CP 16 of reference optical cavities 
|19j for integration times O(10 2 ) seconds as those appropriate for the basic data of the Berlin 
experiment (collected with sets of 10 turntable rotations of 43 seconds each). 

As anticipated, for a further control of the validity of our analysis, we have compared with 
the cryogenic experiment of Ref.|18|. In this case, we have obtained the analogous value 

Q = (13.1 ±2.1) • 1(T 16 Ref.pS] (32) 

from the corresponding Cfc and Sk coefficients. Thus, by using Eq. ([2T)j) and the two values of 
Q reported above, we obtain 

io Xp = (12.0 ± 1.0 ± 3.5) • 1(T 16 Ref.0 (33) 

io Xp = (12.1 ± 1.0 ±2.1) • 10~ 16 Ref.JS] (34) 

where the former uncertainty takes into account the variation of r in Eq. (|28p and the latter 
is both statistical and systematical. 

We emphasize that this stable value of about 12 • 10 -16 is a lower limit for the true 
experimental Aq. At the same time, it could hardly be interpreted as a spurious effect 
of experimental noise. In fact, from the estimates of Ref.|15j. based on the fluctuation- 
dissipation theorem, there is no reason to expect that experiments running in so different 
conditions should exhibit the same experimental noise. 



5. A possible theoretical framework 

Let us now explore a possible theoretical framework that can explain the above experimental 
values. To this end, we shall first use Eq. (|18|) and Eq.Q and re- write the theoretical prediction 
as 

K = - 1)^-/(7, x) (35) 

with 



f(l, x) = y- ~ sin 2 7 cos 2 X~\ cos 2 7 sin 2 xj (36) 

Then, in a flat-space picture of gravity, for an apparatus placed on the Earth's surface the 
vacuum refractive index was estimated to have the value (T| 

2GM n . . 

AA vacuum - 1 ~ -j— ~ 1.4 • 10- 9 37 
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G being Newton's constant and M and R the Earth's mass and radius. Since the meaning of 
"flat-space picture of gravity" might be ambiguous, we shall repeat here the basic derivation. 

The usual interpretation of phenomena in gravitational fields is in terms of a fundamen- 
tally curved space-time. However, some authors |20^ [2T| [22"] have argued that, as light in 
Euclidean space deviates from a straight line in a medium with variable density, an effective 
curvature might also be the consequence of a suitable polarization of the physical flat-space 
vacuum. The substantial phenomenological equivalence with the standard interpretation was 
well summarized by Atkinson as follows [23j : "It is possible, on the one hand, to postu- 
late that the velocity of light is a universal constant, to define natural clocks and measuring 
rods as the standards by which space and time are to be judged and then to discover from 
measurement that space-time is really non-Euclidean. Alternatively, one can define space as 
Euclidean and time as the same everywhere, and discover (from exactly the same measure- 
ments) how the velocity of light and natural clocks, rods and particle inertias really behave in 
the neighborhood of large masses." In this sense, in the flat-space approach, one is adopting 
a "Lorentzian perspective" where physical rods and clocks are held together by the same 
forces underlying the structure of the "ether" (the physical vacuum). Thus the Equivalence 
Principle actually means that the effect of an external gravitational field can be re-absorbed 
into the space-time units of a freely-falling observer so as to preserve local Lorentz covariance. 

Notice that the Equivalence Principle was introduced [23] before General Relativity and, 
as such, does not rely on the notion of a fundamentally curved space-time. Actually, it could 
even be considered misleading within a pure, curved space-time context described by a given 
Riemann tensor [25J. For this reason, it has also been interpreted in terms of an ether flow 
|26j. Nevertheless, regardless of its ultimate physical origin, this fundamental principle implies 
that in a freely-falling frame, given two space-time events that differ by (dx, dy, dz, dt) and 
the local space-time metric 



However, to a closer look, an observer placed on the Earth's surface is equivalent to a 
freely-falling frame up to the presence of the Earth's gravitational field. This leads to the 
weak-field, isotropic modifications of the metric 




(38) 



one gets from ds 2 



the same speed of light as in the absence of any gravitational effect. 



ds 2 = c 2 dt 2 (l 



2GM 



)-(! + 



2GM 



c 2 R 



){dx 2 + dy 2 + dz 2 ) = c 2 dt\ - dl\ 



(39) 
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Here dto and dlo denote respectively the elements of proper time and proper length in terms 
of which, in General Relativity, one would again deduce from ds 2 = the same universal 
value = c. On the other hand, in the flat-space approach the condition ds 2 = is now 
interpreted in terms of a refractive index for the vacuum as in Eq. (|37p and thus light can be 
seen isotropic in only one reference frame, say S. This means that isotropy is only valid if 
the Earth were at rest in the hypothetical S. Otherwise, by introducing the Earth's velocity 
components V£, off-diagonal elements g$i ~ Vi/c in the effective metric are generated by 
the Lorentz transformation for the observer placed on the Earth. These non-zero got can be 
imagined as being due to a directional polarization of the vacuum induced by the now moving 
Earth's gravitational field and express the general property [27] that any metric, locally, can 
be brought into diagonal form by suitable rotations and boosts. By starting in £ with the 
diagonal and isotropic form Eq. ([5U|) . one obtains 

g 0i ~ 2(AA vacuum - 1) — (40) 
c 

and an anisotropy of the speed of light. The value of light anisotropy can also be computed 
by simply composing the isotropic speed c/A/ r va cuum with a given velocity V. In this way, one 
predicts a two-way speed of light as in Eqs.([I]) and (|2]). 



6. Conclusions 

Let us now collect all previous results. For an average colatitude x ~ 38 degrees (as for the 
average location of the Berlin and Dusseldorf laboratories) one finds /(7,x) = (0.60 ± 0.22) 
in the whole range < I7I < it/ 2. Therefore, by using (|37p in Eq. (|35|) . one finds 

Af = (12.6 ± 4.6) • 10~ 16 (V/300 km • s" 1 ) 2 (41) 

in units of the typical speed 300 km/s of most cosmic motions. This estimate is equivalent to 
the previous prediction [28] for the amplitude of the frequency shift measured in a symmetrical 
apparatus (i.e. with two orthogonal rotating resonators) 

^symm = |£| <^> = 2 ^th = ^ g ± q ^ . 1Q -15 ^ 

in terms of the average projected speed yj (v 2 )=(204 ± 36) km/s obtained from a re-analysis 
of the classical experiments. 
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Since these theoretical estimates agree well with the experimental values (f33l) - (l34j) ex- 
tracted from Refs.pJJ [TB], and with the mean amplitude 1.9 • 10~ 15 of the signal measured in 
the symmetric apparatus of Ref.|14j. we conclude that the observed frequency shift, rather 
than being spurious noise of the underlying optical cavities, might also reflect basic properties 
of the vacuum. On the one hand, this appears as a polarizable medium responsible for the 
apparent curvature effects induced by a gravitational field. On the other hand, the observed 
strong random fluctuations of the signal support the view of a stochastic medium, similar to 
a superfluid in a turbulent state of motion. In this sense, the situation might be similar to 
the discovery of the CMBR that, at the beginning, was also interpreted as mere instrumental 
noise. 

A crucial test of our interpretation can be performed in a freely-falling spacecraft. In 
this case, where the vacuum refractive index A/" va cuum for the freely-falling observer is exactly 
unity, the typical instantaneous Ai> should be much smaller (by orders of magnitude) than the 
corresponding 0(1O~ 15 ) value measured with the same interferometer on the Earth's surface. 
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List of Tables 

1 By using Eqs. (J2"3"j) . (J!M|) and (|2"5]). we report the values Qi, their uncertainties 
AQj and the ratio Ri = Qi/AQi for each of the 27 experimental sessions of 
Ref.[3]. These values have been extracted, according to standard error propa- 
gation for a composite observable, from the basic C& and = coefficients 
reported in Fig.2 of Ref.[3j 
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Table 1: By using Eqs. ([23]) . ([24")) and ([25j) . we report the values Qi, their uncertainties AQi 
and the ratio Ri = Qi/AQi for each of the 27 experimental sessions of Ref.[3]. These values 
have been extracted, according to standard error propagation for a composite observable, 
from the basic and Sk = coefficients reported in Fig. 2 of Ref.[3]. 
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